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> ! 1 INTRODUCTION 



The redshifts of lens galaxies in known gravitational lens systems probe the volume 
distribution of lensing mass. Following earlier work by Kochanek, we re-derive the lens 
redshift probability distribution, allowing for mass and number density evolution of 
the lensing galaxies, and apply this test to a much enlarged sample of lens systems. 
From a literature survey of all known lenses, we have selected an unbiased sample of 
15 lenses with complete redshift information. For a flat Universe and no lens evolution, 
we can only put an upper limit on the cosmological constant of £1a < 0.89 at the 95% 
CL. ftA ~ 0.7 and no evolution is consistent with the data. Allowing for evolution in 
an ft m = 0.3, ftA = 0.7 cosmology, we find that the best-fit evolution in c* (i.e., the 
characteristic velocity dispersion in a Schechter-like function) of early- type galaxies, 
in the redshift range z ~ to 1, is rflog [cr*(z)]/dz = —0.10 ± 0.06. This is consistent 
with no evolution and implies that, at 95% CL, of early-type galaxies at z ~ 1 
was at least 63% of its current value. Alternatively, if there is no mass evolution, a 
present-day value of cr* > 175 km s _1 for elliptical galaxies is required (95% CL). 

Key words: cosmology: gravitational lensing - galaxies: general:evolution:mass func- 
tion - quasars: general 



A number of recent measurements all point to a cosmology 
dominated by dark energy. The power spectrum of fluctu- 
ations in the cosmic microwave background indicates a flat 
geometry (e.g., Bennett et al. 2003), while various measure- 
ments of the fraction of the closure density in matter lead 
to consistent values of ft m ~ 0.3 (e.g., Turner 2002). Ob- 
servations of high-redshift supernovae la (Riess et al. 1998; 
Perlmutter et al. 1999) imply that the fraction of the closure 
density in the form of a cosmological constant is ftA « 0.7. 

However, some results obtained from gravitational lens- 
ing statistics appear to be at odds with this picture. The 
incidence of lensing among optically-selected quasars, and 
their image-separation distribution, indicate Q\ < 0.7 at 
95% confidence level (CL; Maoz & Rix 1993; Kochanek 
1996). Similar studies of radio-selected sources (Falco, 
Kochanek, & Munoz 1998) also find, ft A < 0.73 at the 95% 
CL. Chiba & Yoshii (1999) have recalculated the predicted 
lensing statistics, using revised values for the galaxy lumi- 
nosity function parameters, and have argued that a Universe 
with S1a = 0.7 is the most likely one. Kochanek et al. (1998) 
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defended their choice of luminosity function parameters as 
the one that is most consistent with the observed luminosi- 
ties of lens galaxies. Recently, Chae et al. (2002) derived 
constraints on cosmological parameters from analysis of ra- 
dio selected gravitational lenses from the Cosmic Lens All 
Sky Survey (CLASS). They find a best fit ft A = 0.7 ± 0.2 
and constrain ftA > at 95% CL, assuming the most re- 
cent galaxy luminosity function parameters from the 2dF 
survey (Madgwick et al. 2002) and the Sloan Digital Sky 
Survey (SDSS; Blanton et al. 2001). Note that, contrary to 
previous lensing statistics studies, in Chae et al. (2002) the 
velocity dispersion associated with a galaxy of given lumi- 
nosity was left as a free parameter, rather than being set by 
direct observations. Their best fit velocity dispersions have 
relatively low values, which result in low lensing cross sec- 
tions. This may be the main reason why their analysis favors 
a relatively high value of ftA- 

Rix et al. (1994) investigated the impact of hierarchi- 
cal galaxy mergers on the statistics of gravitational lens- 
ing of distant sources. They showed that the total number 
of multiple-image lenses in a sample of potential sources is 
quite insensitive to mergers, but merging leads to a smaller 
mean separation of observed multiple images. Since merging 
does not reduce drastically the expected lensing frequency, 
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it cannot make A-dominated cosmologies compatible with 
the existing lensing observations (e.g., Maoz & Rix 1993). 
Malhotra, Rhoads, & Turner (1997) have suggested that 
small-separation lensing statistics can be reconciled with 
an f^A-dominated Universe by invoking dust in the lensing 
galaxies. The excess number of lensed quasars would then 
be hidden by extinction. Falco et al. (1999) found a mean 
extinction value of AEb-v = 0.04 (0.06) for lens galaxies 
in optically-selected (radio-selected) systems. They tenta- 
tively concluded that the derived total extinction in lenses 
accounts for the small differences in limits on found 
from lensing studies of radio selected vs. optically-selected 
sources. 

Another interesting statistic, that is the focus of this 
paper, was suggested by Kochanek (1992), and is based on 
the redshift of the lensing galaxy in known lens systems. 
Fukugita, Futamase, & Kasai (1990) already showed that the 
expected mean redshift of the deflector increases if 57a dom- 
inates over flrn: due to the large volume at high redshifts. 
Kochanek (1992) obtained a formula for the probability dis- 
tribution of the lens-galaxy redshift given a source (e.g., 
quasar) redshift and the critical radius for lensing. The crit- 
ical radius for lensing was calculated assuming galaxy lenses 
can be represented by singular isothermal spheres (SIS) that 
are isolated (i.e., the lensing is not affected by masses other 
than the lens galaxy) . The Kochanek method is applied by 
comparing for each model (cosmology) the redshift of the 
known lens, zi, with the redshift probability distribution of 
the lens galaxy given the image separation and the source 
redshift, z s . Using a sample of four lenses, Kochanek (1992) 
found that the then-standard cosmologies (without a cos- 
mological constant) are 5-10 times more probable than a 
flat cosmology with Oa 0.9. Since the image separation is 
taken into account as a prior, this method is not affected by 
a bias resulting from the fact that larger separation lenses 
are more easily discovered. This advantage allows, in prin- 
ciple, to use almost all the known lenses, regardless of their 
method of discovery. In §3 we will discuss some exceptions. 

Kochanek (1992) raised several questions and problems 
regarding his method: (i) How to deal with objects in which 
it is difficult to detect the lens directly - should one trust 
a tentative lens redshift based on absorption lines? (ii) Is 
an unevolving galaxy mass distribution justified? (iii) The 
lens system should not be selected based on the lens-galaxy 
redshift. For example, Q2237+0305 was discovered based on 
the quasar emission lines present in the spectra of a nearby 
galaxy. More examples of such objects, which may be found 
in the 2dF survey and the SDSS (e.g., Mortlock & Webster 
2001), arc naturally biased toward lens galaxies at low red- 
shift; (iv) The method requires lenses produced by isolated 
galaxies. Otherwise, the image separations do not necessar- 
ily represent the lens mass, and indeed Kochanek removed 
from his sample objects that are influenced by clusters (e.g., 
Q0957+561). 

The redshift distribution test was re-examined by Hel- 
big & Kayser (1996), who compared the redshifts of six 
lensing galaxies with the probabilities predicted by different 
cosmological models. They assumed that there is no galaxy 
evolution, and that lens galaxies could not be detected be- 
yond a certain magnitude m Jim . For each tested cosmological 
model, they truncated the redshift probability distribution 
beyond the redshift at which the lens galaxy would become 



too faint to have its redshift measured. They concluded that 
the method is not sensitive to cosmological parameters and 
in the foreseeable future it would probably not deliver inter- 
esting constraints. Kochanek (1996) argued, however, that 
the insensitivity found by Helbig & Kayser (1996) was an ar- 
tifact of their statistics, due to the fact that they accounted 
only for lenses with known redshifts, while the fraction of 
the lenses with unmeasurable redshifts was neglected. The 
fraction of lens systems with measured redshifts is smaller 
in a A-dominated Universe. This is due to the higher lens 
redshifts for a given source redshift, combined with a much 
larger luminosity distance to a given redshift. In order to ac- 
count for the selection effect introduced by the detectability 
of the lens galaxy, Kochanek (1996) took into account also 
lens systems for which the lens redshift was not measured. 
Using a total of eight systems (for four of which the lens 
redshift was available) he obtained, for a flat Universe, an 
upper limit of Q.k < 0.9 at the 95% CL. 

Evolution of galaxy mass and number can complicate 
or bias an analysis of the lens redshift distribution. For ex- 
ample, a scarcity of lenses at higher redshifts may be due to 
evolution rather than to the smaller volume of a low Qa cos- 
mology. In fact, there are observational indications for evo- 
lution in the galaxy luminosity function (e.g., Cohen 2002), 
and in the mass-to-light ratio (e.g., Keeton, Kochanek & 
Falco 1998; van Dokkum et al. 2001). Evolution of the lens 
luminosity function and /^-corrections will also affect the 
value of mum- Cohen (2002) studied the galaxy evolution in 
the HDF-North region and concluded that, while the uncer- 
tainties are large, there is no sign of any substantial increase 
in the total mass in stars between z = 1.05 and z = 0.25. 

Keeton (2002) has recently argued that lensing statis- 
tics predictions should be calibrated by counts of distant 
galaxies at z ~ 0.5 — 1 with properties similar to those of 
real lenses, and not by the local density of galaxies. Since 
distant-galaxy number counts depend themselves on cosmo- 
logical volume, when one calibrates the lensing statistics by 
the number counts of distant galaxies, lensing statistics lose 
most of their sensitivity to cosmology. While in principle 
this is true, relating the observed light from high redshift 
galaxies, and the light evolution, to mass is not trivial. The 
advantage of lensing is that it probes mass distribution di- 
rectly. Therefore, it is still a valid approach to compare the 
observed lensing statistics to the predictions of locally cali- 
brated models of lensing populations, and to incorporate the 
mass and density evolution into the models. The mass evo- 
lution indicated by the lensing data can then be compared 
to the luminosity evolution implied by number counts, to 
gain insight both on cosmology and on the relation between 
mass and light. 

The purpose of this paper is to revisit the question of 
the lens redshift distribution, its dependence on cosmology, 
and on the evolution of the lens population, for the larger 
sample of lensed quasars now available. In §2 we rederive 
Kochanek's formalism for the expected lens redshift distri- 
bution, but allow for number and mass evolution of the lens 
population. We also investigate the properties of the lens 
redshift distribution, and its sensitivity to different param- 
eters. In order to define an unbiased sample, in §3 we con- 
duct a literature survey for known gravitational lens system 
and list them (in an Appendix) with their basic parameters. 
In §4 we constrain the cosmological and mass-evolution pa- 
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rameter space. We discuss the results in §5, and give a short 
summary in §6. 



2 CALCULATION OF THE LENS REDSHIFT 
DISTRIBUTION 

2.1 Formalism 

The differential optical depth to lensing per unit redshift is 

where n(9, z) is the comoving number density of lenses that 
have critical angular radius for lensing 9 at redshift z, S 
is the cross section for lensing, and cdt/dz is the proper 
distance interval. 

We assume that the lens galaxies can be represented by 
a SIS potential. The SIS assumption is consistent with lens 
data, galaxy dynamics, and the X-ray emission from ellip- 
ticals (e.g., Fabbiano 1989; Rix et al. 1997; Treu & Koop- 
mans 2002). Moreover, it is known that ellipticity in the lens 
galaxy mainly affects the relative numbers of two- and four- 
image lenses but not their overall measures (e.g., Keeton, 
Kochanek, & Seljak 1997). The critical radius (Einstein ra- 
dius), and cross section, for lensing in a SIS potential are 
given by 

2 Di. - 
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where Di, D 3 , and Di s are the angular diameter distances 
between observer-lens, observer-source, and lens-source re- 
spectively (e.g., Fukugita et al. 1992). Note that in the SIS 
model, the separation between lensed images is always 29. 
Unless mentioned otherwise, we use "filled-beam" angular- 
diameter distances throughout this paper (as opposed to 
"empty-beam", [Dyer & Roeder 1972, 1973]). Je is a free 
parameter that represents the relation between dark halos 
and their stellar component, as explained below. 

The velocity dispersion odm of the mass distribution 
and the observed stellar velocity dispersion a need not be 
the same. Turner, Ostriker, & Gott (1984) argued that if el- 
liptical galaxies have dark massive halos with more extended 
distributions than those of the visible stars, this dark mate- 
rial must necessarily have a greater velocity dispersion than 
the visible stars, and a factor of (3/2) is required to give 
an r~ 2 density distribution (cf. Gott 1977). Kochanek (1992) 
adopted a parameter Je that relates the dark-matter veloc- 
ity dispersion odm and the stellar velocity dispersion a: 



0~DM = /eCT. 



(4) 



Dynamical models by Franx (1993), and Kochanek (1993, 
1994) demonstrated that the assumptions leading to the 
(3/2) 1/2 factor were incorrect. However, White & Davis 
(1996) argued based on X-ray observations that there is a 
strong indication that dark matter halos are dynamically 
hotter than the luminous stars. In our work, we have kept 
/_b as a free parameter, since it mimics the effects of: (i) 
systematic errors in the velocity-dispersion measurements 
of galaxies; (ii) large scale structures that change the typ- 
ical separation between images (e.g., Martel, Premadi, & 



Matzner 2002; see §3 below) - an increase of an arbitrary or- 
der F in the velocity dispersion is equivalent to an increase 
of F 2 in the typical separation 9 (i.e., 9 oc <r 2 ); (iii) softened 
isothermal sphere potentials which tend to decrease the typ- 
ical image separations (e.g., Narayan & Bartelmann 1996), 
and could be represented by fs somewhat smaller than 1. 
Kochanek (1996) already noted that small core radii will 
lead to only slight modifications of the tail of the lens red- 
shift probability distribution. Moreover, there is no evidence 
for cores in lensing galaxies (e.g., Rusin & Ma 2001). 

Martel et al. (2002) compared the distribution of image 
separations obtained from ray-tracing simulations in CDM 
models (Premadi et al. 2001), with analytical calculations. 
They found that the presence of the background matter 
tends to increase the image separations produced by lens- 
ing galaxies. However, they found this effect to be small, of 
order of 20% or less, and independent of the cosmological pa- 
rameters (see also Bernstein & Fischer 1999). Bromley et al. 
(1998; see also Christlein 2000, Zabludoff & Mulchaey 1998) 
showed that the luminosity function of early type galaxies 
depends on environment and that richer environments tend 
to have a higher ratio of dwarf to giant galaxies than the 
field. Keeton, Christlein, & Zabludoff (2000) showed that 
this effect nearly (but not entirely) cancels the effect of the 
background matter, making the distribution of image sep- 
arations essentially independent of environment. They pre- 
dicted that lenses in groups have a mean image separation 
that is ~ 0.2" smaller than that of lenses in the field. It is 
possible that all these factors can affect the images separa- 
tion by up to ±20%. This can be mimicked by introducing 
(0.8) 1/2 < f E < (1.2) 1 / 2 . 

In order to derive the lens redshift probability we 
need the velocity-dispersion distribution function of galax- 
ies. Such a function is not yet available directly, although it 
is expected that the SDSS and 2dF surveys, combined with 
a good understanding of aperture effects, will produce one 
soon. Instead, we start with the Schechter (1976) function, 



n(L, z) — n*(z) 



L*{z) 



exp - 



dL 



L*(z) J L*(z) ' 



(5) 



which provides a good description of the galaxy luminosity 
distribution. Using the SDSS data, Blanton et al. (2001) find 
that this function fits the number density of SDSS galaxies 
better than an alternative non-parametric model. In order 
to relate the luminosity to velocity dispersion, we use the 
Faber- Jackson relation (Faber & Jackson 1976; Kormendy & 
Djorgovski 1989; Bernardi et al. 2001) for early-type galax- 



(6) 



L,(z) \a*(z) 



A similar relation (i.e., a Tully-Fisher relation), but with 
different cr„ and 7, is adopted for spiral galaxies. 
From Equations 6 and 2 we obtain 



7 /2 



dL 7 



L,(z) \9,(z) J ' L*(z) 2 \9*(z) 
where 

a 1 \ A ( o~*{z)\ 2 Di s 2 
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In our rederivation of the lens redshift probability, we 
allow for evolution in n», L», and cr,, with the following 
parametrization. 



n„(z) 


= n,10 Pz , 


L*{z) 


= LA0 Qz 


a.(z) 


= a*W Uz , 



where P, Q, and U are constants. U represent an evolution 
in the characteristic cr* , of a Schechter-like "mass" function. 
Note that the definitions of P and Q are somewhat different 
than the ones used by Lin et al (1999). The definition of Q 
is the same as the one used by Cohen (2002). Throughout 
the paper we use n„, L», cr*, and #» to refer to the values at 
zero redshift, as opposed to n, (z), £»(z), &*{z), and 0* (z). 
Our parametrization is consistent with hierarchical galaxy 
merging. Mergers change only the position of the galaxies 
on the Faber- Jackson relation but not the basic shape of 
the relation. The scatter in the Faber- Jackson relation (e.g., 
Kormendy & Djorgovski 1989; Bernardi et al. 2001) can be 
taken into account using a Monte-Carlo simulation, as we 
discuss in §2.3. 

Combining Equations 7 and 5, to obtain n(6,z), and 
inserting n(8,z), along with the evolution parametrization 
(Eq. 9, 10, and 11) into Eq. 1, we obtain the optical depth 
per unit redshift for a system with image separation 26 and 
source redshift z s , 

x(l + ,)»^?£(£)** 1+ " )+1 (12) 
xexp[-(£)* 7 10-«^], 
where the normalization, 

r N =4ttV| (^) 2 . (13) 

The lens light emission (as opposed to the lens mass) does 
not affect lensing, and therefore Q does not appear in Equa- 
tion 12. For a given lens system, the dependence of dr/dz on 
z gives the relative probability of finding the lens at different 
z's. 

2.2 Choices for the Lens Parameters 

Throughout this paper we adopt the Schechter function pa- 
rameters for the fej-band, found by Madgwick et al. (2002) 
from a sample of 75, 000 galaxies in the 2dF survey with 
median redshift of z ~ 0.1 (see also Blanton et al. 2001 for 
a similar result based on a smaller sample using the SDSS). 
We assume a — —0.54 ±0.01 for early-type galaxies (ellipti- 
cals and SO) and a = — 1.16±0.01 as the average of late-type 
galaxies (weighted by their abundance). The relative abun- 
dance of early-type to late-type galaxies is adopted from the 
local density, as obtained by Madgwick et al. (2002), while 
the relative abundance of ellipticals (E) to SO is obtained 
from Postman & Geller (1984; E:S0=39:61). We thus adopt 
n, = 1.46 x 10~ 2 ; 0.61 x 10~ 2 ; 0.39 x 10~ 2 h 3 Mpc" 3 for 
spirals, SOs, and ellipticals, respectively. The typical uncer- 
tainties in these number are of order 5%. Note that these 
values are different from those used by Kochanek (1992, 
1996). We do not incorporate evolution in a, as the red- 
shift distribution test is not very sensitive to a (see §2.3). 







Spiral 


SO 


Elliptical 


a 




-1.16 


-0.54 


-0.54 


ra* 


[h 3 Mpc- 3 ] 


1.46 x 1(T 2 


0.61 x 10~ 2 


0.39 x 10~ 2 


7 




2.6 


4.0 


4.0 


cr * 


[km s -1 ] 


144 


206 


225 



Table 1. The default parameters for the luminosity function and 
the Faber-Jackson relation. 



A constant a is also justified by observations (Cohen 2002). 
For the Faber-Jackson relations, Bernardi et al. (2001) find 
7 = 4.00, 3.91, 3.95, 3.92 in the SDSS g, r, i, and z-bands, 
respectively, with about 5% error. The spectral dependence 
of the slope 7, is smaller than the typical scatter in the 
Faber-Jackson relation (~ 40% in a ln-normal distribution; 
Bernardi et al. 2001). We therefore adopt 7 = 4.0 for early- 
type galaxies. For the Tully-Fisher relation of spirals, we 
assume 7 = 2.6, following Kochanek (1992). We take cr* 
from Fukugita & Turner (1991), based on de Vaucouleurs 
& Olson (1982), and Efstathiou, Ellis, & Peterson (1988): 
cr, = 225 km s _1 for ellipticals; cr* = 206 km s _1 for SOs; 
and <7, = 144 km s _1 for spirals. Note that M* (the mag- 
nitude of a cr, galaxy) given by Madgwick et al. (2001) and 
that given by Efstathiou et al. (1988) agree well, after apply- 
ing the transformation between the b.j and Bt magnitude 
systems (i.e., bj = B T + 0.29). 

Table 1 lists the default parameters we will assume in 
this paper. Unless mentioned otherwise, fs = 1, fi m = 0.3, 
Q A = 0.7, U = 0, P = 0. 

2.3 What Can Be Tested by the Lens Redshift 
Distribution 

In this section we investigate the (generally complicated) 
behavior of dr /dz{9, z 3 ) (Equation 12) as a function of the 
parameters involved in the problem. We also study the ef- 
fect of the scatter in the Faber-Jackson relation on the lens 
redshift distribution. 

Kochanek (1992) noted that the lens redshift distribu- 
tion cuts off more sharply at high redshift than the total 
optical depth dr/dz (i.e., Eq 12 integrated over all f?'s), be- 
cause the critical radius constraint introduces an exponen- 
tial term from the Schechter function beyond the redshift 
at which an L, galaxy is required to produce the observed 
separation. The sharp cutoff makes dr / dz{9, z s ) a powerful 
tool. For example, the existence of relatively high redshift 
lenses, even in a small sample, can reject models in which 
such redshifts are beyond the cutoff in the lens-redshift dis- 
tribution. 

Figure 1 shows the lens redshift probability distribu- 
tion for several different values of a given parameter, while 
keeping the rest of the parameters fixed. Unless otherwise 
noted, we use the morphological mix of galaxies with the 
default parameters described in §2.2, along with z s — 2 and 
26 — 1". Since it is not practical to scan all the parameter 
space shown in Figure 1, this figure can help us to isolate 
the parameters that have the most profound effect on the 
lens redshift distribution. 

Evolution Parameters U and P: As seen in Fig 1-g 
the parameter U, which describes the evolution of the char- 
acteristic galaxy velocity dispersion, has a considerable ef- 
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Figure 1. The lens redshift probability (normalized to one) density distribution as a function of redshift for various choices of parameters. 
The vertical axis is the probability in arbitrary units. In all panels, unless otherwise noted, we use the default parameters: z s = 2, 29 = 1", 
n m = 0.3, Qjy = 0.7, Je = 1, U = 0, P = 0, and the morphological mix with the parameters given in Table 1. 



feet. Assuming P ~ 0, negative U decreases the most prob- 
able redshift of the lens z distribution, and makes the distri- 
bution narrower. This suggests that redshift evolution can 
impact the results considerably, and can therefore be con- 
strained by the data. The effect of the number density evo- 
lution parameter, P (see Fig 1-h), is somewhat smaller, and 
mainly changes the most probable redshift, as opposed to U 
that affects also the cut-off redshift at which the lens-redshift 
probability goes to zero. 

Cosmology: From Figure 1-d, it appears that the lens 
redshift distribution test in a flat Universe is sensitive to Ha 
mainly when Qa <; 0.8. 

Morphological type: It is well known (e.g., Maoz & 
Rix 1993) that almost all the contribution to lensing opti- 
cal depth comes from early-type galaxies. Since the lensing 
galaxy type is not well constrained, in the analysis we will 
take it to be a mix of all the galaxy types ' . As seen in Fig- 



t Although several lens galaxies are clearly disk systems (e.g., 
B 1600+434), even in those cases it is not clear whether they are 
SO's or later types. 



ure 1-e, as long as early-type galaxies dominate the lensing 
galaxy population, there will be no dramatic effect on the 
final results. 

fE: As discussed in §2, the /e parameter (Figure 1- 
c) can mimic the uncertainties in the relation between im- 
age separation and mass, or a systematic error in a,. Al- 
though the effect of complex environments on the image 
separation is not clear enough, it seems that in most cases 
it cannot change the image separation by more than 20% 
(see discussion in §1). Therefore, for reasonable values of 
(0.8) 1/2 < f E < (T2) 1/2 , f E can somewhat affect the results 
of this work, but not as considerably as the mass evolution 
or a large cosmological constant. However, a systematic er- 
ror of more than about 10% in u* can have some effect on 
the lens redshift distribution. 

Angular diameter distance: We have tested the ef- 
fect of the angular diameter distance formulae (filled-beam 
and empty-beam approximations) on the lens redshift distri- 
bution (not shown in Fig 1). Using empty-beam approxima- 
tion somewhat decreases the typical lens redshift relative to 
the filled-beam approximation, but again, the effect is con- 
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Figure 2. The lens-redshift probability distribution for the de- 
fault parameters, and z s = 2 and 26 = 1" (solid line). The dottcd- 
dashed and dashed lines were calculated using the same param- 
eters, but for the mean of 10,000 Monte-Carlo realizations with 
In- normal scatter introduced to cr/cr*, of 20% and 40% respec- 
tively. 



siderably smaller than that of mass evolution or of a large 
cosmological constant. 

7: Even a 10% change in 7 does not cause significant 
variations in the lens redshift distribution. Note that the 
effect of 7 in Figure 1-i is tested only for U = 0, and its 
effect on the exponential term of Eq. 12 is therefore limited. 

a: The effect of changing a is shown in Figure 1-f. Madg- 
wick et al. (2002) give a 2% error on a. Changing a by 10% 
(not shown) has a negligible effect, even smaller then the ef- 
fect of changing 7. Setting a to a completely different value 
(e.g., a = —1.1; Efstathiou et al. 1988), somewhat shifts the 
peak of the lens-redshift probability distribution. 

An important issue that needs to be considered is the 
scatter in the Faber- Jackson relation, which could alter the 
results. The observed scatter in the Faber- Jackson relation 
is about 40% in a ln-normal distribution (Bernardi et al. 
2001). To study the effect of such scatter, we have per- 
formed Monte-Carlo simulations. In each realization we per- 
turbed cr/a* by multiplying it by a factor randomly selected 
from a ln-normal distribution with a given standard devi- 
ation. We performed two sets of Monte-Carlo simulations 
with standard deviation of 20%, and 40%. For each set we 
calculated the lens-redshift probability, averaged over the 
realizations, as a function of redshift. Figure 2 shows the 
mean effect of the scatter in the Faber- Jackson relation on 
the lens-redshift probability distribution. The solid line is 
the unscattered lens-redshift probability distribution as a 
function of redshift for the default parameters. The dashed- 
dotted line is the same, but when 20% scatter is introduced 
to a /it*, and the dashed line is for a 40% scatter. It is clear 
that the Faber- Jackson scatter makes the distribution some- 
what broader. However, it does not change the distribution 
significantly. Note that using the lens redshift distribution 
function based on a velocity dispersion distribution function 
would overcome this problem. 

Given the most probable values of the parameters, the 



mass evolution U has the most important effect on the lens- 
redshift distribution. This is not surprising since U appears 
in the exponential term of Equation 12. Next in importance 
are the cosmological constant, the number density evolution, 
P, and /e- In what follows, we therefore investigate the 
parameters found to be important to this problem, namely 
U, P, f E , and Qa- 



3 SAMPLE SELECTION 

Using the lens redshift distribution of known lenses to con- 
strain cosmological parameters and galaxy evolution re- 
quires a sample free of selection effects and biases. In Ta- 
ble Al, we have compiled a list of 71 galaxy lens systems 
we are aware of (July 2002) with their basic properties. 
The main source for this list is the CASTLES database^ 
(Munoz et al. 1999) that, at the time of writing, included 
65 systems. We have added six additional systems discov- 
ered recently: FIRST J100424.9+122922 (Lacy et al. 2002); 
HS0810+2554 (Reimers et al. 2002); PMNJ0134-0931 (Winn 
2002b); FIRST J0816+5003 (Lehar et al. 2001); HE0435- 
1223 (Wisotzki et al. 2002); PSS2322+1944 (Carilli et al. 
2002). 

There are several lens-system candidates that, although 
listed in the CASTLES database, are of uncertain nature, 
usually due to differences in the spectra of the supposed 
lensed images or the non-detection of a lensing galaxy. Re- 
jection of these candidates based on the absence of a lens- 
ing galaxy is problematic. For example, a lensing galaxy 
in a real lens system may be undetected because it is dis- 
tant and hence faint. Systematically rejecting such systems 
would bias the sample against detecting the effect of a cos- 
mological constant. We therefore reject lens candidates only 
if their image spectra show pronounced differences. Based 
on this criterion, Q0252-3249 (Morgan et al. 2000) is prob- 
ably not a real lens. In any event, this object has z s = 2.24, 
and we will see below that even if it were a real lens system 
it would not enter our sample. 

As mentioned above, some lens systems have been dis- 
covered based on the properties of the lens galaxy, thus pre- 
ferring low-redshift lenses. On the other hand, it is diffi- 
cult to detect lensed quasars for which the lens galaxy is 
nearby, based on quasar selection. After removing from our 
sample lens systems discovered based on galaxy properties, 
we therefore expect a bias against low- 2 lenses. In princi- 
ple, we should correct for this bias by modifying the lens- 
redshift distribution. In practice, the probability for such 
lenses at 2 <J 0.1 is very low anyway and we can safely ig- 
nore this correction. Lenses rejected based on their discovery 
method (i.e., lens-galaxy selected, as opposed to source se- 
lected) are: Q2237+0305 (Huchra et al. 1985); CFRS03.1077 
(Crampton et al. 2002); FIRST J0816+5003 (Lehar et al. 
2001). We have also rejected all the lens candidates dis- 
covered in the HST medium deep survey (Ratnatunga, 
Griffiths, & Ostrander 1999): HST14176+5226; HST12531- 
2914; HST14164+5215; HST15433+5352; HST01247+0532; 
HST01248+0531; HST16302+8230; HST16309+5215; 
HST12369+6212; HST18078+4600; and in the Hubble Deep 

t http://cfa-www.harvard.edu/castles/index.html 
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Fields: 

HDFS2232509-603243 (Barkana, Blandford, & Hogg 1999). 
These lenses were selected based on their morphology (e.g., 
the detection of a possible lensing galaxy) and therefore, 
could introduce a bias. Rejecting all these lens systems does 
not introduces any bias, since we are effectively ignoring 
entire surveys. Furthermore, only two of these systems have 
complete redshift information, and as we will see below, their 
source redshift is too high to be included in our sample. In 
Column 10 of Table Al these objects are marked "L". 

There are several lens candidates (listed as binary 
quasars in the CASTLES database) whose nature is not 
clear (e.g., Q1634+267, Peng et al. 1999; Q2345+007, Small, 
Sargent, & Steidel 1997, Green et al. 2002). Even if these sys- 
tems are real lenses, we reject them from our sample, since 
their separation is too large for galaxy lensing, and they 
could be affected by cluster-mass objects. 

Kochanek (1996) accounted for systems with missing 
redshifts by considering the probability that one cannot de- 
tect the lens galaxy given the limiting magnitude to which 
a system has been probed. However, the lens galaxy bright- 
ness depends on galaxy luminosity evolution and this could 
complicate the analysis. In order to avoid this problem we 
have defined a sample that is source-redshift limited. We 
have selected all systems with z s < 2.1 and image sepa- 
rations smaller than 4" (Sample I). This z s selection crite- 
rion has several advantages: (i) Even for a de-Sitter Universe 
(Qa = 1), the most probable redshift is < 1.2 and therefore 
we can probe the evolution factors within a similar redshift 
range to that studied by optical surveys (e.g., Cohen 2002); 
(ii) For z 3 < 2.1 we can expect that the redshifts for most of 
the lenses have been measured, since they are not too high. 
However, selecting a limit (as we have) that is based on the 
largest z s for which z t is measured for most systems, is some- 
what biased. For example, in a A-dominated Universe, the 
typical zi will be higher and therefore the maximum z s (for 
which the sample is complete) will be lower. To check for 
the robustness of the final results, we have defined a second 
sample (Sample II) for which the maximum z s is 1.7. 

The lens redshift information in the samples is not 
complete. There are three lens systems for which the lens 
redshift is not available yet. These are J0158-4325 (Mor- 
gan et al. 1999), HE0435-1223 (Wisotzki et al. 2002), and 
HS0810+2554 (Reimers et al. 2002). All three systems were 
discovered recently, and spectral observations have not been 
attempted. However, the lens galaxies are detectable and are 
actually quite bright (R = 19.4, i — 18.1, and / ~ 19, re- 
spectively) relative to other lens galaxies. Therefore, drop- 
ping these systems from our sample will not introduce a bias. 
For all but one of the lensing galaxies in samples I and II, 
the lens redshift has been measured spectroscopically. The 
exception is FBQ0951+2635, for which we have used the 
BVRI photometry of Schechter et al. (1998), corrected for 
Galactic extinction (Schlegel, Finkbeiner, & Davis 1998), to 
calculate a photometric redshift by \ minimization relative 
to galaxy spectral templates from Kinney et al. (1996). We 
find a photometric redshift of 0.23, 0.24, and 0.28, for Sb, 
SO and E0 galaxy types, respectively. Kochanek et al. (2000) 
showed that early- type field galaxies with < z < 1, as rep- 
resented by lens galaxies in lensed quasar systems, lie on 
the same fundamental plane as those in rich clusters at sim- 
ilar redshifts. They used this to estimate the "fundamental- 



plane redshift" of lens galaxies by matching the lens galaxy 
effective radius, magnitude and velocity dispersion to the 
fundamental plane of rich clusters. The Kochanek et al. 
(2000) fundamental-plane redshift for FBQ0951+2635, al- 
though somewhat dependent on cosmology through the lu- 
minosity distance, is consistent (zi = 0.21 ± 0.03) with our 
photometric redshift. Based on these estimates, we have 
adopted, zi = 0.25, for FBQ0951+2635. In §4 we check for 
the robustness of this assumption by changing the redshift 
of this object by ±0.05. The lens system B2045+265 was ex- 
cluded from all the samples, since its source redshift is most 
likely in error (see Fassnacht et al. 1999). 

The image separation criterion of 4" was introduced in 
order to remove lenses that are influenced by complex en- 
vironments (e.g., clusters). The separation is used as prior 
information in calculating the lens redshift probability for 
each system. Therefore, even if we reject a lens that has a 
large separation, only because its lens redshift is very low 
relative to the source redshift, this selection will not bias 
the experiment. On this basis we have removed Q0957+561, 
and RXJ0921+4529 from Samples I and II. To test for the 
robustness of this cut, we have defined Sample III, in which 
we have used sample I plus these two ~ 6" separation lenses. 
Samples I, II, and III contain 15, 11, and 17 sources, respec- 
tively. 

For six lenses in sample I for which an elliptical isother- 
mal sphere plus shear model is available in the CASTLES 
database, the critical radius (9) from the model is about half 
the observed image separation (or half the maximum sep- 
aration in quadruple systems) to 10% accuracy or better. 
Therefore, we adopt half the image separation as represent- 
ing the critical radius in all systems. 



4 ANALYSIS 

Figure 3 shows the lens redshift probability distribution, as 
calculated using Equation 12 with the default parameters, 
for each of the lenses in samples I, II, and III. In each panel 
the name of the lens along with its source redshift (z s ), lens 
redshift (z;), and image separation (29) in arcseconds are 
listed. The actual lens redshift is marked by an arrow on 
the bottom axis. 

Our statistical analysis is based on the maximum like- 
lihood technique. For each lens system (i), given the model 
(described by Equation 12), we calculate the probability 
density Pi(zi\z 3 , 9; X), which is normalized to one, where 
X are the model parameters (e.g., 51a, U), and z s and 9 are 
priors. The log likelihood (In C) as a function of X is given 
by 

ln£ = ^lnP;. (14) 

i 

Testing one parameter at a time (i.e., Xj), the confidence 
interval (CI) within a given confidence limit (CL), is the 
region found below the peak of the log likelihood by 0.5, 
2.0, and 3.3174, for 68%, 95%, and 99% CLs, respectively 
(e.g., Press et al. 1992). 

Figure 4 shows ln£ for the full sample I. The dotted 
horizontal lines are the 68%, 95%, and 99% CLs, from top 
to bottom, respectively. From this figure we see that the 
most likely cosmological model (assuming no evolution) is 
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Figure 3. Lens redshift probability distributions calculated for morphological mix of galaxies with the default parameters (see §2.2), 
for each of the 17 lenses in samples I, II, and III. The vertical axis is the probability in arbitrary units. In each panel the name of the 
lens along with its source redshift (z s ), lens redshift (z;) and image separation in arcseconds (28) are noted. The observed lens redshift 
is marked by an arrow on the bottom axis. 



Q A = 0.36±£:?o> with upper limits Q A < 0.89 at the 95% 
CL, and fl A < 0.95 at the 99% CL. 

The maximum-likelihood technique is not free of bias 
(i.e., the most likely value could be biased; e.g., Efron 1982; 
Shao & Tu 1995). To test the results for bias and possible 
outliers, we use the Quenouille-Tukey jackknife technique 
(e.g., Efron 1982; Efron & Tibshirani 1993). For a sample 
of iV "observations", we remove one "observation" (i.e., one 
system) at a time, and calculate In £ as a function of Q,a, 
for the N — 1 remaining observations. The jackknife esti- 
mate for the bias, B, in the parameter r\ (i.e., Qa) is given 
by B = (N — 1)(< r\i > — ij). Where 7/ is estimated using N 
observations, and r/i is the same, but using N — 1 observa- 
tions at a time, where i is the index of observation dropped. 
This bias should be subtracted from r\ to obtain the cor- 
rected r). We find a bias of about A^a ~ —0.2. Taking into 
account the bias, the most probable £l A « 0.6 is well within 
the la error and therefore we will ignore the bias. 

Sample II favors a somewhat larger Q A = 0.601°;^. 



The higher Q A favored by sample II may be due to the z s 
cut we introduced (see §3) . However, the Q A values from the 
two samples are consistent within the uncertainty indicated 
by the likelihood ratios. Sample III favors an even higher 
cosmological constant, Q A — 0.77^q 22 (1°" errors). This re- 
sult is not surprising, since cluster-affected lenses which have 
large separation will tend to have lenses at redshifts higher 
than expected, assuming low Q, A and ignoring the effect of 
the cluster. This will increase the derived value of the cos- 
mological constant. 

As noted in §3, the redshift of the FBQ0951+2635 lens 
is based on a photometric redshift. In order to test for the 
robustness of this datum, we have recalculated In C when 
perturbing the lens redshift within the plausible range of 
redshifts (zi = 0.25 ± 0.05; see §3), and have found that it 
has little effect. 

Given a model, the maximum-likelihood technique gives 
information about the best fit parameters of the model. 
However, the technique does not give information on the ac- 
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Figure 4. Likelihood plot for Sample I. The solid curve is the log 
of the likelihood as a function of Qa for + £l\ = 1- The rest 
of the parameters arc the default ones. The dotted lines are the 
la, 2a and 99% CLs. The vertical lines on the bottom axis mark 
the points where In C crosses the dotted lines. 
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Figure 5. Likelihood contours in the U-P plane. The solid curves 
are the, Icr, 2cr and 99% CL lines for sample I, where the plus sign 
mark the position of maximum likelihood. The dotted curves are 
these contours for sample II, for which the triangle marks the 
position of the maximum. 
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tual validity ("goodness of fit") of any model. To this end, 
we have performed a set of Monte-Carlo simulations. Given 
a set of model parameters, in each simulation we used the 
source redshift (z s ) and image separation (29) of the objects 
in sample I to generate a set of 1000 "sample I" realizations. 
Each realization contains 15 lens redshifts (like sample I). 
We then calculated the likelihood expectation value, and its 
distribution, as a function of Oa- We find that the likelihood 
expectation value is about 4.4 ± 1.1 (la errors), with some 
dependence on S7a, in good agreement with the maximum 
likelihood shown in Figure 4. We conclude that the lenses in 
sample I have redshifts consistent with the probability func- 
tion given by Equation 12, and that our basic assumptions 
are probably valid. 

Next, we set the cosmology to fi m = 0.3, S7a = 0.7, 
as found by Bennett et al. (2003), and maximize the likeli- 
hood for the evolution parameters. Note that the evolution 
parameters U and P are applied to all morphological types. 
As the cross section for spiral lensing is considerably smaller 
than for ellipticals, U and P effectively trace the evolution in 
early- type galaxies. Figure 5 shows the likelihood contours 
in the U-P plane. The solid curves are the la, 2a and 99% 
CL contours (for two degree of freedom, see Press et al. 1992) 
for sample I, where the plus sign marks the position of max- 
imum likelihood. The dotted curves mark the same contours 
for sample II, for which the triangle marks the position of 
the maximum likelihood. We can see that both samples give 
consistent results. As expected, the lens redshift distribu- 
tion is a better probe of U than of P. Marginalizing over P, 
Figure 6 shows ln£ as obtained by integrating C(P, U) over 
P (i.e., In J C(P,U)dP). In the absence of prior knowledge 

regarding P, we obtain U = -OAOto^to^-oll^ where 
the errors are for the la, 2a and 99%-CL CIs, respectively. 
Similarly, marginalizing over U, we get P = +0.71^2 (not 
shown). Finally, Figure 7 shows, for Q. m — 0.3, S1a = 0.7, 
and no evolution in the number density of lenses (P = 0), 
the likelihood contours in the U-/e plane. Assuming that 



6 




Lens a. Evolution U 

Figure 6. Marginalized likelihood as obtained by integrating 
C(P,U) over P (i.e., In J C(P,U)dP). Symbols as in Fig 4. 

the dark-matter velocity dispersion and the stellar velocity 
dispersion are the same, and that large scale structures do 
not change considerably the image separations, then Je rep- 
resent a constant factor that multiplies the <j,'s in Table 1. 
Thus, it mimics a systematic error in the "canonical" value 
of it* used in this paper. Assuming no evolution (P = 0, 
U = 0), we obtain a best fit f E = 0.90in?',l u '.i2 to.is ( lcr > 
2a and 99%-CL errors). 



5 DISCUSSION 

The lens redshift distribution test, like other cosmological 
tests, depends on assumptions, some of which may not be 
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Figure 7. Likelihood contours in the U-Je plane. The solid 
curves are the, la, 2a and 99% CL lines for sample I, where 
the plus sign mark the position of maximum likelihood. 



correct. Based on the analysis we have conducted in §2.3, we 
would say that our most troubling assumption is that im- 
age separation truly represents the galaxy mass. Large scale 
structure can influence the image separation, and therefore 
we are probing some complex combination of galaxy mass 
and intervening structures and not the mass of the galaxy on 
its own. This is important as we use the "mass function" of 
galaxies to derive the lens redshift distribution. Our current 
knowledge suggests that in most cases, the image separation 
is not affected by more than 20% (see §2.1) by large scale 
structure. 

Extinction can influence somewhat the number statis- 
tics of lenses, by lowering the observed number of optically 
selected lenses. Falco et al. (1999) found a mean value of 
AEb-v = 0.05 for lensed galaxies. The lens redshift distri- 
bution test is not influenced by dust, due to the fact that 
known lens systems are used and their image separations are 
taken as priors. However, for calculating the lens-redshift 
probability, we are using the observed mix of galaxy mor- 
phological types. Since spirals are more dusty, there could be 
a selection-induced lower fraction of spiral lenses relative to 
elliptical lenses, and this effect would lower the cr*-evolution 
parameter, U, (or Qa) deduced by the lens redshift distribu- 
tion method. However, the lens population is dominated by 
early- type galaxies by virtue of their larger mass. As shown 
in Figure 1, lowering the fraction of spiral lenses will not 
affect the predictions dramatically. 

We showed that the lens redshift distribution test is 
sensitive to high values of the cosmological constant (e.g., 
SIa ^ 0.8). Using the lens redshift distribution, and assum- 
ing no evolution (U — 0, P = 0) and a flat Universe, the 
current sample puts an upper limit of Q.\ < 0.89 at the 95% 
CL, and Q A < 0.95 at the 99% CL. This is similar to the up- 
per limit obtained by Kochanek (1996) for a much smaller 
sample, using the same method. Indeed, we find that the 
observations are consistent with the "concordance" value of 
0,a ~ 0.7. Note, however, that this method with the current 
sample is still not a sensitive probe of Qa for fl A <J 0.8; our 
results are equally consistent with Qa = 0. 



We have shown that the lens redshift distribution is a 
sensitive probe of the velocity-dispersion evolution of early- 
type galaxies. This arises from the exponential dependence 
on the velocity-dispersion evolution parameter U in Equa- 
tion 12. Assuming £l m = 0.3, 57a = 0.7 and no prior knowl- 
edge of U or P, then the best fit values are: P — +0.7t\\t 
and U — —0.10 ± 0.06. Our results apply to the redshift 
range ~ 0.1 to ~ 0.9. The lower limit is set by the median 
redshift of galaxies used to obtain the 2dF luminosity func- 
tion (Madgwick et al. 2002) and the Faber-Jackson relation 
(Bernardi et al. 2001). The upper limit is set by the highest 
redshift lens in sample I (z = 0.93). Our results are con- 
sistent with no evolution. Positive U (mass increases with 
redshift) can probably be excluded as physically improb- 
able. Our measurements are then useful for setting a lower 
limit on this "mass" evolution parameter. Negative U would 
mean that the typical characteristic velocity dispersion (in 
a Schechter-like function) a,, was lower in the past. Such 
evolution is expected in the hierarchical galaxy merging sce- 
nario, as the typical galaxy was less massive in the past. 
Our results imply that, at 95% CL, <r* of early-type galaxies 
at z ~ 1 was at least 63% of its current value. As we have 
shown, the lens redshift distribution does not strongly con- 
strain the number- density evolution P. If we take as a prior 
P — —0.26, as obtained by Im et al. (2002) and supported 
by Schade et al. (1999) from counts of early- type galaxies to 
z « 1, we see from Figure 5 that U = -OMto^t^to^l, 
where the errors are the la, 2a and 99% CL errors, respec- 
tively. This places a more severe limit on the allowed evolu- 
tion in cr»; at 95% CL, a* at 2 fa 1 was at least 77% of its 
current value. 

Most studies of galaxy evolution attempt to constrain 
the evolution of the light, or of the mass-to- light ratio de- 
fined by (M/L)(z) = (M/L) W Tz , rather than the mass 
evolution. The luminosity evolution, number evolution, and 
mass-to-light ratio evolution have been studied in various 
works that produced several estimates for the evolution pa- 
rameters P, Q, and T. For example, Cohen (2002) studied 
the luminosity function of galaxies in the HDF-North region, 
in the redshift range to 1.5. She found that galaxies with 
strong absorption lines become brighter with 2 (Q ~ 0.6 in 
all bands) , while galaxies with detectable emission lines show 
a smaller change in L„, with Q ~ 0.3. She also found that 
the comoving number density and stellar mass in galaxies is 
approximately constant out to 2 ~ 1, and with more uncer- 
tainty, to z ~ 1.3. Keeton, Kochanek, & Falco (1998) com- 
bined the photometry and lens modeling of 17 lens galaxies 
in the redshift range 0.1 to ~ 1. They found that the lens 
galaxy B-band mass-to-light ratio decreases with redshift 
with T « -0.5 ± 0.1 (for Q. m = 0.1). In Table 2 we list 
these, and several other evolutionary studies of early-type 
galaxies, and their best fit parameters. The similar values 
of — T and Q found by most studies, around 0.6, suggest 
considerable luminosity evolution with little mass evolution, 
consistent with our findings. Note that the evolution param- 
eters listed in Table 2 are determined for early-type galaxies 
in different environments (field and clusters). Interestingly, 
the evolution parameters are similar. 

In a recent paper, Davis, Huterer & Krauss (2002) use 
the CLASS and Jodrell-Very Large Array Astrometric Sur- 
vey lensing statistics to constrain the characteristic velocity 
dispersion of elliptical galaxies. This is done using the joint 
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Reference 


Zmax 


Parameters 


van Dokkum et al. (1998) 


0.83 


T « 


-0.45 


Keeton, Kochanek, & Falco (1998) 


~ 1 


T as 


-0.5 


Lin et al. (1999) 


0.55 


P = 


-0.7 ± 0.2 






Q = 


0.8 ±0.2 


van Dokkum et al. (2001) 


0.55 


T = 


-0.59 ± 0.15 


Cohen (2002) 


~ 1 


Q~ 


0.6 


Rusin et al. (2002b) 


~ 1 


T = 


-0.56 ± 0.04 


Im et al. (2002) 


~ 1 


Q = 


0.8 ±0.3 






p = 


-0.26 ±0.20 



Table 2. Evolutionary parameters found using other methods. 
Zmax is the maximum redshift for which the parameters were 
fitted. Since the parameters are cosmology dependent, we choose 
only works which assume a "standard" cosmology (i.e., Q m = 0.3, 
Qa = 0.7), or cosmologies that have luminosity distances that 
are not too different (e.g., 50 = 0.1 as used by Lin et al. 1999, or 
Qq = 0.1 used by Keeton et al. 1998). Evolution in luminosity, 
number, and mass-to-light ratio are parameterized by Q, P, and 
T, respectively. 



likelihood for the number of lenses, the redshift distribu- 
tion of lensing galaxies, and the angular-separation distri- 
bution of lensed images. Their analysis gives 138 km s' 1 < 
cr* < 206 km s _1 (95% CL). However, their lens-redshift 
distribution analysis include only the six out of 12 lenses 
in their sample for which redshift information is available. 
For a given image separation, lenses with higher redshift 
have higher masses. Therefore, failing to take into account 
the missing redshifts will bias cr, toward lower values. Fur- 
thermore, the "canonical" value of cr, (e.g., 225 km s _1 for 
ellipticals; Fukugita & Turner 1991) is measured at redshift 
zero, while the lensing analysis constrains cr* at the mean 
redshift of the lenses. Therefore, it is hard to tell if the anal- 
ysis by Chae et al. (2002) and Davis et al. (2002) is detecting 
evolution in cr* or an error in the old value of cr* . The mean 
lens-galaxy redshift of sample I is z7 = 0.54. Our best fit 
value of U = —0.10 can be mimicked by a constant cr* that 
is about 88% (= W Uz ') of the local cr* we have used in this 
work. However, we can see from Figure 7 that a cr* that is 
<; 80% of the values we have used in this paper (see Table 1), 
can be rejected since it would require positive evolution (U) 
at the 95% CL. Assuming no evolution (U = 0, P = 0), we 
have found that the most probable f E = 0.90lo;o?'+o.i2,i"o'.i5 
(la, 2a and 99%-CL errors). Therefore, a present-day value 
of a, > 175 km s _1 for elliptical galaxies is required (95% 
CL). This value is in agreement with a constant velocity dis- 
persion as found by Chae et al. (2002), cr* = 198±|? km s" 1 
for ellipticals. Note that Chae et al. (2002) use a = —1.0, as 
opposed to a = —0.54. 

Thus, a constant value of a, « 200 km s _1 or a mildly 
evolving value that has approximately this average would 
explain most of the current lensing statistics data in an Q.\ = 
0.7 cosmology. Since lensing cross section is proportional to 
cr 4 , such a value of cr, may lower the lensing cross section 
in lens surveys by a large enough amount to explain the 
small number of observed lensed systems despite the large 
volume associated with a cosmological constant. At the same 
time our results show that such cr, is still large enough to 
reproduce the observed redshift distribution of known lenses. 



6 SUMMARY 

We have compared the observed and predicted redshift dis- 
tributions of lensing galaxies. Expanding upon the work of 
Kochanek (1992, 1996) we have used a much enlarged obser- 
vational sample and incorporated possible mass and number 
density evolution of lens galaxies. The lens-redshift distribu- 
tion in a flat Universe is sensitive to the large volume out 
to a given redshift in the presence of a large cosmological 
constant (i.e., Qa ii 0.8), and especially to evolution of the 
characteristic velocity dispersion of the lensing population. 
The sensitivity of this test to velocity-dispersion evolution is 
unique, as it probes the total mass of stars and dark matter 
within ~ 10 kpc, without any assumptions on a possible bias 
factor between mass to light. We have explored the robust- 
ness of our inferences to various types of deviations from our 
basic assumptions, such as errors in the luminosity-function 
parameters, the lens model, and the scatter in the Faber- 
Jackson relation. Unless there are gross errors in our knowl- 
edge of these parameters, they have a small effect on our 
results. 

There are now 71 lenses known in the literature and we 
have compiled their basic properties (Appendix A). From 
these, we have selected a sample of 15 lenses with separa- 
tions smaller than 4" and source redshift z 3 < 2.1, along 
with two other test samples. Assuming no mass evolution 
(0 & z ^ 1) of early-type galaxies, we can put an upper 
limit of D.A < 0.89 at the 95% CL. On the other hand, as- 
suming an Q, m — 0.3, Qa = 0.7 cosmology we find that the 
best fit characteristic velocity dispersion (a,) and number 
evolution of early- type galaxies, in the range z ~ to 1, 
are dlog 10 a*(z)/dz = —0.10 ± 0.06 and dlog 10 n*(z)/dz = 
+0.7tl;2. Our find ings are consistent with no evolution in 
the early-type galaxy population and limit the allowed in- 
crease of the characteristic cr, between z ~ 1 and the 
present. Finally, we have considered the influence of chang- 
ing the estimate of the present-day value of cr* . We find that 
a velocity dispersion that is <; 80% of the standard value of 
cr* would require positive velocity-dispersion evolution (i.e., 
galaxies were more massive in the past) at the 95% CL, and 
therefore can be rejected. A value of a, as 200 km s _1 for 
early- types, with no evolution out to z « 1, is consistent 
with our data. Such a value of cr, could possibly also re- 
solve the discrepancy between the results of previous lensing 
statistics studies and an Qa = 0.7 cosmology. 
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Name 


R.A. 


Dec. 


Zs 




separation 






Grade 


Sample 


Ref 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 




(8) 


(9) 


(10) 
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HDFS2232509-603243 
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: 22 
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22,27 


PSS2322+1944 
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4.12 




1.5 
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61 



Table Al. The columns arc: (1) lens name; (2) R.A. [h:m:s], J2000.0; (3) Dec. [d:m:s], J2000.0; (4) source redshift, z s ; (5) lens redshift, z !; 
(6) image separation in arcseconds, 28; (7) lens magnitude; (8) Number of images, R for ring; (9) Grade for the likelihood that the object is 
a lens, as adopted from the CASTLES database: A=I'd bet my life, B=I'd bet your life, and C=I'd bet your life and you should worry; (10) 
sample: L-lens was discovered based on the lens-properties (rather than the source), I- Sample I in this paper, II- Sample II in this paper, 
Ill-additional objects to those in Sample I, that are found in sample III; (11) References. 
List of references: 

1 - Lacy ct al. (2002); 2 - Warren ct al. (1996); 3 - Ratnatunga, Griffiths, & Ostrander (1999); 4 - Phillips et al. (2000); 5 - Morgan ct al. 
(1999); 6 - Wisotzki ct al. (1999); 7 - Reimers et al. (2002); 8 - Morgan et al. (2000); 9 - Crampton et al. (2002); 10 - Falco, Lchar, & Shapiro 
(1997); 11 - Gregg et al. (2000); 12 - Fassnacht & Cohen (1998); 13 - Marlow ct al. (2001); 14 - Tonry & Kochanck (1999); 15 - Ibata ct al. 

(1999) ; 16 - Kobayashi et al. (2002); 17 - Koopmans et al. (1999); 18 - Winn et al. (2002a); 19 - Tonry & Kochanek (2000); 20 - Winn et al. 

(2000) ; 21 - Winn ct al. (2001); 22 - Rusin ct al. (2001a); 23 - Barkana, Blandford, & Hogg (1999); 24 - Winn et al. (2002b); 25 - Gregg et al. 
(2002); 26 - Siemiginowska ct al. (1998); 27 - Lubin ct al. (2000); 28 - Lchar et al. (2001); 29 - Augusto et al. (2001); 30 - Kncib, Cohen, & 
Hjorth (2000); 31 - Burud ct al. (2002b); 32 - Kochanek et al. (2000); 33 - Morgan et al. (2001); 34 - Burud ct al. (2002a); 35 - Murioz et al. 

(2001) ; 36 - Rusin ct al. (2002a); 37 - Surdcj ct al. (1997); 38 - Marlow et al. (1999); 39 - Fassnacht et al. (1996); 40 - Rusin ct al. (2001b); 41 
- Sykcs ct al. (1998); 42 - Koopmans & Treu (2002); 43 - Wisotzki et al. (2002) ; 44 - Fassnacht et al. (1999); 45 - Myers et al. (1999); 46 - 
Lidman et al. (2000); 47 - Lehar ct al. (1993); 48 - Hall et al. (2002); 49 - Huchra ct al. (1985); 50 - Langston et al. (1989); 51 - O'Dea et al. 
(1992); 52 - Wiklind & Combes (1995); 53 - Schechter et al. (1998); 54 - Lehar ct al. (2000); 55 - Walsh, Carswell, & Weymann (1979); 56 - 
Hewett et al. (1994); 57 - Eisenhardt ct al. (1996); 58 - Tonry (1998); 59 - Magain ct al. (1988); 60 - Wiklind & Combes (1996); 61 - Carilli et 
al. (2002); 62 - Cohen, Lawrence, & Blandford (2002); 63 - Hagen & Reimers (2000) 

Notes: 

a - Based on absorption line spectrum (Gregg et al. 2000). 

b - Associated with a group at z = 0.29 and virial mass of ~ 300 km s _1 (Fassnacht & Lubin 2002). 
c - Associated with a small group and possibly a second group at z = 0.52 (Tonry & Kochanek 1999). 
d - Two late-type lensing galaxies (Koopmans et al. 1999). 

e - Redshift is based on the Kochanek et al. (2000) method in which one requires that the galaxy properties be consistent with the passively 

evolving fundamental plane of early- type galaxies (Winn et al. 2002a). 

/ - Evidence for a foreground group of galaxies at z = 0.34 (Tonry & Kochanck 2000). 

g - Redshift is based on the Kochanek et al. (2000) fundamental-plane method, plus the possible identification of [OIII] in emission (Winn et 
al. 2002b). 

h - The most probable redshift is of an Mg II absorber at z = 1.1349 (Siemiginowska et al. 1998). 
i - Complex lens, probably by two galaxies at z = 0.3157 and z = 0.5883 (Augusto et al. 2001). 
j - Massive cluster at z = 0.769 with a ~ 840 km s" 1 (Kncib, Cohen, & Hjorth 2000). 

k - Spectral observations reveal a rich metal-line absorption system consisting of a strong Mg II doublet and associated Fc I and Fc II absorption 
features, all at an intervening redshift of z = 0.684, suggestive of a lensing galaxy (Morgan ct al. 2001). 

I - By cross correlating the spectrum with galaxy templates, Burud et al. (2002a) obtain a tentative redshift estimate of z = 0.50, in agreement 
with the fundamental-plane method estimate (Kochanck ct al. 2000). 

m - Possibly a member of an X-ray cluster (z = 0.32) centered on the field (Murioz et al. 2001). 

n - Based on photometric redshift (this paper). The fundamental-plane method (Kochanck et al. 2000) gives z; = 0.21 ± 0.03. 
o - z; = 0.41 ± 0.05 is based on Kochanek ct al. (2000) fundamental-plane method. 
p - Based on lensing galaxy colors z = 0.4 ±0.1 (Wisotzki ct al. 2002). 
q - Cluster at z = 0.36 (e.g., Chae 1999). 

r - Based on the probable identification of the lens galaxy 4000Abreak in the quasar spectra (Lacy et al. 2002). 
s - Based on the fundamental-plane method (Kochanek et al. 2000) z = 0.88 ± 0.07. 

t - Absorption lines at z = 2.344 and Mgll absorption at z = 1.085 (Surdej et al. 1997). Kochanek et al. (2000) based on the fundamental-plane 
method give z = 0.78 ± 0.07. Therefore, we have adopted z = 1.085. 

u - Marlow et al. (1999) argue that the lens galaxy is sub-L„, so it has z > 0.5, or is highly reddened. 
v - Lensing by a compact group of three z f» 1 galaxies (Rusin et al. 2001b). 
w - Possible cluster at z ~ 1.7 (Kncib et al. 1998). 

x - Keeton & Kochanek (1997) find that the lens time-delay cannot be fit by an isolated lens galaxy, but that it can be well fitted by including 
a perturbation from the nearby group of galaxies. Tonry (1998) measures the group velocity dispersion of 330 km s . 

y - Eisenhardt et al. (1996) estimate the lens-galaxy redshift as z = 0.9 ± 0.3. Goodrich ct al. (1996) find strong evidence for Mgll absorption 

at z = 1.316 in the spectrum of the QSO and weaker evidence for a possible continuum and absorption system at z = 0.893. Lacy, Rawlings, 

& Serjeant (1998) identify an Ho emission line with z = 0.914. The Kochanek et al. (2000) fundamental plane method gives z = 0.78 ± 0.05. 

z - Based on the Kochanck et al. (2000) fundamental plane method. 

aa - L. Koopmans and T. Treu, private communication (radio ring detected in radio). 

ab - Excess of galaxies at z ~ 0.9 (Faurc et al. 2002). 
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